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1 Introduction 


This is a final report for NASA grant NAG2-513, “The analysis of control 
trajectories using symbolic and database computing.” This report covers 
the period January 1, 1992 through December 31, 1994. This report also 
comprises the formal semi-annual status reports for this grant for the pe- 
riods June 30-December 31, 1993, January 1 — June 30, 1994, an d June 
1 — December 31, 1994- 

The research supported by this grant is broadly concerned with the 
symbolic computation, mixed numeric-symbolic computation, and database 
computation of trajectories of dynamical systems, especially control sys- 
tems. The NASA Technical Officer for this grant is Dr. George Meyer, 
NASA Ames Research Center, Mail Stop 210-3, Moffet Field, California, 
94035. In Section 2, we provide an overview of the progress made during 
the report period. Section 3 contains bibliographic references for the articles 
related to this grant that were completed during this period. 

2 Review of Work During Report Period 

Trajectories and Approximating Series. Nonlinear dynamical systems, con- 
trol systems, and hybrid systems all have trajectories which describe the 
behavior of the system as time evolves. The analysis of these trajecto- 
ries and the derivations of specialized numerical algorithms which approx- 
imate them is greatly aided by the exact computation of symbolic series 
which approximate the trajectories. These are series involving differential 
operators — developing the algorithmic foundations of the symbolic analysis 
of such series has been an important focus for our work. Taylor series, Lie 
series, and Fliess-Chen series are all examples of such series. 

Technically, it has turned out to be useful to start with a ring of functions 
R (such as polynomials or formal power series, which can be managed easily 
symbolically) and to view differential operators as higher order derivations 
of R. A derivation is simply a linear map of R to itself satisfying Leibnitz’ 
rule, while a higher order derivation can be viewed for the purposes here as 
a formal symbolic expression in derivations of R. 

The Cayley algebra of trees. Symbolic computations involving differen- 
tial operators typically involve cancellations. It turns out that developing 
efficient algorithms for computing symbolically with differential operators 
depends upon exploiting these cancellations [Crouch92]. In prior work, we 
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developed a data structure involving trees labeled with derivations to keep 
track of certain of these cancellations. By doing this, one can develop more 
efficient algorithms for the symbolic computation of differential operators. 
The analysis of these cancellations is simplified by imposing an addition and 
multiplication on trees so that the action of the differential operators fac- 
tors through an action of the algebra of trees. We call the resulting algebra 
the Cayley algebra of trees, which we denote by T. It turns out that it is 
a Hopf algebra. These and related ideas are described in [Grossman92b], 
where a number of algorithms are presented for the symbolic manipulation 
of expressions involving derivations. 

Actions of differential operators. An application in which the symbolic 
computation of derivations is important is the derivation of numerical al- 
gorithms, especially application-specific numerical algorithms. One means 
of deriving numerical algorithms in an invariant or basis-free fashion is to 
compute the action of a flow (which can be considered as a differential op- 
erator) on a test or observation function. The symbolic computation of this 
action requires extending the ideas above involving trees from the context 
of algebras to the context of modules. More precisely, we showed how to 
view the ring of observation functions R as a T-module. These ideas are 
described in [Crouch91] and [Crouch93]. In [Crouch92], we specialize this 
calculus to systems with polynomial coefficients and show how it can be 
used to derive first integrals for such systems. In [Doffou95], we specialize 
this calculus to derive differential invariants for vector field systems with 
polynomial coefficients. This provides an explicit symbolic algorithm for a 
problem which has been outstanding for many years. In [Grossman94b] we 
applied closely related ideas to study normal forms of control systems. 

Geometrically stable integration algorithms. State spaces G of nonlin- 
ear systems often satisfy geometric constraints: in general, numerical algo- 
rithms do not preserve these constraints. We used the Cayley algebra of 
trees to derive specialized numerical algorithms with the property that by 
construction they preserve certain geometric structures. It turns out that 
one of the key technical ideas is to work with the ring R of all observations 
R = {/ : G — ► R} of the state space G rather than with the state space G 
itself. This approach allows numerical algorithms to be derived in an invari- 
ant or basis independent fashion [Crouch93]. Specific numerical algorithms 
were derived by using the Cayley algebra of trees and viewing R as a T- 
module. In particular, we derived application- specific numerical algorithms 
for a rotating rigid body (G = 50(3)) and a ball rolling on an inclined plane 
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(<j = E( 3)) in [Crouch95]. 

Hybrid Systems. Discrete systems, such as automata, and continuous sys- 
tems, such as control systems, have both internal representations involving 
states, and external representations involving inputs and outputs. The ex- 
ternal representation can be represented by a formal series or generating 
series p. Not all such series are associated with systems. A realization the- 
orem gives conditions on p which describes the system externally so that 
there' exists an internal representation involving states. 

We have already mentioned how Hopf algebras arise and can be used 
to organize combinatorial computations, especially those that arise in sym- 
bolic computation. In [Grossman92a], we show how there is a natural Hopf 
algebra H associated to a nonlinear system with inputs and outputs and 
how the formal series p can be viewed els element of the dual H* . The Hopf 
algebra H acts on the space of all observations R of the system. 

This led to a new class of realization theorems. In the case of nonlinear 
systems with input and outputs the realization theorem becomes the Fliess 
Theorem [Grossman92a]; in the case of finite automata, the realization the- 
orem becomes the Myhill-Nerode theorem [Grossman95a]. 

This suggested that hybrid systems be defined through an internal de- 
scription involving observation functions R rather than through an internal 
description involving states. From this viewpoint [Grossman95a], a hybrid 
system consists of a commutative algebra R and a Hopf algebra II , together 
with an action of H on R. The action expresses the dynamics. The extreme 
cases of a nonlinear control system and an automaton are recovered by as- 
suming special forms for H . A hybrid system simply consists of a general H, 
which corresponds to a collection of interacting nonlinear control systems 
and automata. In [Grossman95c] these ideas were used to reveal a structural 
similarity between hybrid systems and quantum automata. 

Trajectory Stores. An interesting class of hybrid systems arises by consid- 
ering a collection of control systems, each one corresponding to a separate 
mode of the hybrid system. Transitions between modes are specified by 
an automaton. These are mode-switched hybrid systems. Each state of a 
finite automaton corresponds to a nonlinear control system describing an 
underlying physical system in a particular regime. Such hybrid systems are 
a special case of the definition above. An important open problem for such 
systems is: Given a sequence of points in state space, find a path connecting 
them using flows of the underlying control systems and satisfying the mode 
changes allowed by the underlying automaton. This is the path planning 
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problem. 

Trying to compute such paths in real time is usually prohibitively ex- 
pensive, except in special cases. On the other hand, one can store many 
precomputed trajectory segments and try to select an appropriate sequence 
of them which approximate the desired path connecting the sequence of 
given points. This is the approach we took [Grossman93b]. It involves 
managing large collections of trajectory segments or TrajectoryStores using 
object managers. 

PTool. Our solution to path planning for hybrid systems is a variant of the 
classical approach of trading space for time. In our case, we traded high per- 
formance computing for high performance data management. In particular, 
over several years, we developed a persistent object manager called PTool 
[Grossman94a], [Grossman95b], and [Grossman95d] which supports the type 
of high performance data management required in scientific computing. Un- 
like, a traditional database management system, PTool is optimized for the 
type of data access for applications such as the TrajectoryStore described 
above. Recently, we have developed variants of PTool suited for clustered 
computing environments. In [Grossman94c], we described how PTool could 
be used for hybrid systems involving rules and agents. 

Other Activities. We organized a workshop “Application- Specific Symbolic 
Techniques in High Performance Computing Environments” at the Fields 
Institute in Waterloo, Ontario from October 18-20, 1993. The workshop was 
funded by $7500 provided by NASA grant NAG 2-513 and 87500 provided 
by the Fields Institute. 

Approximately fifty scientists, engineers, and computer scientists at- 
tended the workshop, which focused on three themes: 

1. using symbolic methods to derive application specific numerical algo- 
rithms 

2. using symbolic methods to analyze nonlinear systems 

3. algorithms and software for hybrid systems. 

3 Publications Written During the Report Period 

[Crouch92] P.E. Crouch and R. L. Grossman, “The Explicit Computation 
of Integration Algorithms and First Integrals for Ordinary Differential 
Equations With Polynomial Coefficients Using Trees,” Proceedings of 
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the 1992 International Symposium on Symbolic and Algebraic Compu- 
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